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摘 x. 我 国 西北 干旱 区 内 陆 河 流域 水 资源 匮乏 ,水 资源 利用 主要 用 于 农业 生产 ,准确 估算 内 陆 河 流 
域 蒸 散 发 与 农业 灌溉 效率 ,对 研究 内 陆 河 流域 气候 变化 和 水 资源 合理 利用 具有 重要 作用 。 利 用 基于 
地 表 能 量 平 衡 方程 的 SEBAL 模 型 ,对 2017 一 2018 年 芷 勒 河流 域 蒸 散 量 进行 定量 估算 与 时 空 分 布 特 
征 分 析 , 并 结合 降水 量 与 净 灌 溉 水 量 数据 ,对 朴 勒 河流 域 昌 马 灌区 的 年 内 灌溉 水 有 效 利 用 系数 进行 
估算 。 结 果 表 明 :(1) 聋 勒 河流 域 2017 一 2018 年 日 均 ET 呈 单 峰 变 化 趋势 ,最 大 值 为 6 月 的 5.03 mm: 
d !, 最 小 值 为 12 月 的 0.55 mm: d ,并 存在 明显 的 空间 分 布 差 异 。(2) AARRE ETARDA, 
夏季 ET 达到 最 高 的 201. 83 mm, 春 秋 次 之 ,冬季 最 低 为 53. 92 mm;ET 由 东南 向 西北 逐渐 减 小 ,流域 
上 游 E7 明显 高 于 中 下 游 地 区 。(3) 昌 马 灌区 各 灌溉 时 段 不 同 的 灌溉 水 量 造成 了 各 灌 季 蒸 散 量 的 差 
异 ,灌区 有 7 高 值 区 主要 分 布 在 中 部 与 东南 部 , 低 值 区 主要 分 布 在 西北 部 和 灌区 边缘 。(4) 昌 马 灌区 年 
内 灌溉 水 有 效 利 用 系数 呈 下 降 趋 势 ,其 中 春灌 、 夏 灌 、 秋 灌 和 冬 灌 分 别 为 0.76.0.71.0.69 和 0. 55, = 


均 灌 流水 有 效 利用 系数 为 0.67。 


X Ga. 蒸 散发 ; SEBAL 模型 ， 疏 勒 河流 域 ; 灌溉 水 有 效 利用 系数 ; 昌 马 灌区 


文章 编号 : 


蒸 散发 (87) 是 流域 水 文人 循环 和 水 量 平衡 的 主 
要 环节 ,作为 区 域 能 量 平衡 的 重要 组 成 部 分 , 它 反 
上 映 了 在 大 气 .土壤 和 植被 变化 的 条 件 下 ,水 和 能 量 
通 量 之 间 的 相互 作用 "””。 因 此 ,通过 葵 散 发 的 研 
究 可 深入 了 解 水 资源 的 形成 过 程 与 变化 规律 ,对 水 
资源 评价 与 干旱 监测 具有 重要 作用 *“ ,尤其 是 在 
西北 干旱 区 , 蔡 散 发 的 准确 估算 和 监测 对 于 区 域 水 
资源 管理 与 利用 意义 重大 。 

FAROE A FE SEM Ti EA ARB SOE BEL — 
能 量 平 衡 法 s BEA ASG ATA RE GE SE? S, 
但 均 存 在 工作 量 大 ,数据 获取 的 有 限 性 以 及 不 能 
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以 合理 的 成 本 对 大 尺度 区 域 进行 全 面 观测 等 局 
限 性 。 不 断 发 展 的 遥感 技术 手段 为 定量 估算 大 
斥 度 区 域 的 地 表 蒸 散 量 提供 了 方法 ,结合 遥感 数 
据 和 模型 对 区 域 蒸 散发 的 反 演 与 应 用 也 成 为 国 
内 外 研究 的 重要 方向 。 目 前 广泛 使 用 的 计算 模 
型 有 陆 面 地 表 能 量 平衡 方程 (SEBAL) 、 地 表 能 量 
平衡 系统 (SEBS) .基于 地 表 温 度 和 植被 指数 的 梯 
形 空间 模型 . 双 源 能 量 平衡 模型 (TSEB) .METRIC 
方法 等 "-。 其 中 ,SEBAL 模 型 具有 明确 的 物理 基 
础 ,具有 适应 性 广 精度 高 的 特点 ,可 用 于 大 尺度 区 
域 的 长 时 间 地 表 蒸 散 量 计算 ,因而 在 国内 外 得 到 了 
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广泛 的 应 用 ,如 国外 学 者 RAHIMI SANTOS 
ELNMER 等 "利用 不 同 数 据 源 结合 SEBAL 模 型 进 
行 蒸 散 发 反 演 ,并 在 流域 尺度 上 取得 较 好 的 反 演 结 
果 ;国内 学 者 李 宝 富 呈 、 刘 春雨 "9 VS EMT. 
CHANG 等 “分别 以 中 国 西北 干旱 区 为 研究 区 , 利 
用 SEBAL 模 型 进行 蒸 散 发 反 演 研 究 ,结果 表明 SE- 
BAL 模 型 能 够 有 效 实现 对 干旱 区 ET IB. AL, 
本 文 将 采用 SEBAL 模 型 探究 玖 勒 河流 域 的 蒸 散 发 
特征 。 

内 外 众多 学 者 对 灌溉 效率 的 评价 指标 与 内 
涵 开 展 了 很 多 研究 ,针对 灌溉 水 有 效 利 用 系数 也 做 
了 相应 的 探讨 。 借 助 遥 感 手段 可 快速 获得 灌区 
蒸 散发 ,将 灌溉 地 消耗 的 水 量 作为 灌溉 水 有 效 利 用 
量 并 用 于 灌溉 效率 的 估算 ,避免 了 传统 灌溉 效率 测 
算 的 复杂 性 和 存在 的 监测 困难 ,为 农业 水 资源 高 效 
利用 提供 了 技术 支撑 ””。 国 内 外 学 者 利用 遥感 落 
散发 模型 实现 了 对 灌区 的 耗 水 分 析 与 灌溉 评价 , 结 
果 表 明基 于 遥感 蒸 散 发 的 灌溉 水 利用 效率 评价 方 
法 可 快速 准确 的 实现 灌溉 效率 评价 ,具有 较 强 的 可 
操作 性 2 >), LATO PLA AK ARES. ,主要 存在 
有 限 的 水 资源 与 急剧 增加 的 人 口 .土地 以 及 与 流域 
内 农业 生态 用 水 的 矛盾 3, 使 得 流域 内 水 资源 的 
合理 配置 显得 尤为 重要 。 农 业 用 水 约 占 琉 勒 河流 
域 每 年 水 资源 量 的 89% ,用 水 方式 主要 依赖 于 水 利 
工程 的 灌溉 ™, 对 流域 内 农业 用 水 效率 的 快速 准确 
评价 与 提高 也 显得 尤为 迫切 。 因 此 明确 流域 实际 
蒸 散 发 的 时 空格 局 和 灌溉 用 水 效率 ,对 于 该 区 水 资 
源 可 持续 利用 具有 重要 作用 。 

目前 对 下 勒 河流 域 蒸 散 发 研究 主要 为 利用 传统 
方法 的 测定 和 基于 遥感 模型 的 短 时 间 尺 度 反 演 呈 -加 ， 
XT eat AD Ta THE X OE BE Se ZE EE TB hl EA BAR 
等 方面 ” ,人 鲜 有 对 年 内 蒸 散发 的 变化 特征 以 及 
对 灌区 灌溉 效率 在 不 同 灌溉 季节 的 分 析 评 价 研 
究 。 本 文 利 用 MODIS 数据 和 气象 观测 数据 ,采用 
SEBAL 模型 定量 估算 了 踊 勒 河流 域 2017 一 2018 
年 内 的 地 表 营 散发 ,明确 琉 蒜 河流 域 年 度 ET AE 
化 特征 。 结 合流 域 昌 马 灌区 各 灌溉 季节 的 灌溉 
水 量 数据 ,对 昌 马 灌区 年 内 灌溉 水 有 效 利 用 系数 
进行 了 估算 ,以 便 更 为 真实 地 反映 区 域内 实际 的 
灌溉 效率 ,为 玻 勒 河流 域 水 资源 管理 和 分 配 等 提 
供 参 考 。 


1 数据 与 方法 


1.1 研究 区 概况 

玻 勒 河流 域 地 处 我 国 西北 干旱 区 , 介 于 38?00' 
~ 42°48' N,92°11' ~ 98°30’ 下 之 间 。 流 域 面积 约 为 
1. 25x10? km ,是 河西 走廊 三 大 内 陆 河流 域 之 一 。 
臣 勒 河流 域 上 游 为 高 海拔 山区 ,发 育 有 典型 大 陆 性 
冰川 ,并 分 布 有 多 年 冻 土 ; 流 域 中 下 游 地 势 低 平 , 荒 
漠 分 布 广泛 ,土地 沙化 严重 “。 玖 勒 河流 域内 土地 
利用 类 型 主要 为 低 覆 盖 度 草地 `. 沙 地 和 车 漠 ,植被 
分 布 稀少 ,生态 环境 脆弱 。 玻 勤 河 流域 属 典型 的 干 
旱 荒 漠 性 气候 ,具有 降水 稀少 蒸发 强烈 的 气候 特 
点 。 玖 勒 河流 域 年 均 降雨 量 低 于 60 mm, 多 年 平 
均 气 温 6.9 ~ 8.8 °C ,年 蒸发 量 最 高 可 达 3 000 mm, 
属于 严重 干旱 地 区 。 玻 勒 河流 域内 辖 昌 马 、 双 
塔 . 花 海 三 大 灌区 , 玻 勒 河 灌区 历经 多 年 的 水 利 建 
设 与 生态 保护 综合 规划 项 目的 实施 ,灌溉 效率 显 
著 提 升 , 并 被 指定 为 甘肃 省 灌溉 用 水 有 效 利用 的 
FE b BE DX 
1.2 数据 来 源 

研究 所 用 数据 主要 有 遥感 数据 .DEM 数 据 AC 
象 数 据 、 土 地 利用 覆 被 数据 和 灌溉 统计 数据 。ET 反 
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Fig. 1 Land use of Shule River Basin and the location of 


Changma irrigated area 
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演 模 型 输入 的 各 地 表 参 数 由 美国 NASA 网 站 
(https ://ladsweb. nascom. nasa. gov/search/) 中 的 MO- 
DIS 系列 产品 获得 。 本 文 根 据 晴 天 条 件 和 云 量 、 影 
像 质 量 选择 玻 勒 河流 域 2017 年 3 月 ~2018 年 2 月 对 
应 的 MODIS 遥感 影像 。 气 象 数据 来 源 于 中 国 气 象 
科学 数据 共享 服务 网 (http ://ede. ema gov. cn/home. 
do) ,选用 卫星 过 境 时 对 应 地 面 观测 数据 ,包括 气温 、 
风速 .相对 湿度 、 气 压 等 地 表 常 规 气象 数据 。DEM 
数据 为 空间 分 辨 率 为 30 m 的 ASTER GDEM V2 产 
品 ,来 源 于 地 理 空间 数据 云 (http://www. gscloud. 
cn/)。 土 地 利用 枚 被 数据 来 源 于 中 国 科学 院 资 源 环 
境 科学 数据 中 心 (http://www. resde. cn/)。 用 于 计算 
TE BUSCAR BO ER DC TRE CIBUS OR UST HEN E i TT Bat 
域 水 资源 管理 局 (http ://slt. gansu. gov. cn/slhglj/) 中 
提供 的 灌溉 管理 数据 。 
1.3 SEBAL 模型 计算 原理 

SEBAL 模 型 是 基于 遥感 的 陆 面 能 量 平 衡 模型 ， 
通过 遥感 数据 与 气象 参数 对 各 能 量 分 量 进行 估算 ， 
得 到 区 域 的 实际 蒸 散发 ”。 区 域 的 陆 面 能 量 平衡 
方程 为 : 


AET-R,-G-H (1) 

IUP: AET 是 潜 热 通 量 (Wm 7); R, 是 净 辐 射 通 
量 (W*m 7); G 是 土壤 热 通 量 (Wem); H 是 感 热 
iat (W*m 7); 和 是 潜 热 蒸发 系数 (W.m .mm ')。 
1.3.1 45:548 ER, 

R,-(1-a)R, * R, - R,-(1-2)R, (2) 
SUH: a 是 地 表 反 照 率 ; R 是 太阳 短波 总 辐射 ; UR, 
是 大 气 长 波 辐射 ; R, 是 地 面 长 波 辐射 ; e 是 地 表 比 
辐射 率 。 
1.3.2 土壤 热 通 量 G 
a x (0.003 8a+ 0.007 4a?) us 


G- / 
a h - 0.98 x NDVI‘) 


WP: 7, 为 地 表 温 度 ; NOVI 为 归 一 化 植被 指数 。 
1.3.3 ARMS 

_ paso dT 
Poh 
RP: pw 是 空气 密度 (kg.m-?); C, 是 空气 定 压 
比热容 ( 取 1004J.kg-…K- ;ad7 是 不 同 高 度 之 
间 的 温度 差 ; ru, 是 空气 动力 学 阻抗 (sm ')。 其 
中 mr， 和 dT 通过 选取 冷 热 像 元 进行 循环 欠 代 计 
算 取 得 。 


H (4) 


. 


1.3.4 i& Xv 58iEdK X HT ILU 
的 为 地 面 瞬时 数据 ,通过 蒸发 比 不 变 的 方法 将 瞬时 
蒸 散 量 扩展 到 日 蒸 散 量 汪 ,计算 公式 如 下 : 
__AET (5) 
RC 


ET. = R 24 X €x 864 00 

* [2.501 - 0.002 361 x (T, - 273.15)] x 10° 
A: o 为 蒸发 比 ; ETa AA ARH; Ra NH 
Seti at (Wm 7); T, 为 地 表 温 度 (K)。 
1.4 蒸 散 发 时 间 尺 度 扩展 及 验证 
1.4.1 参考 作物 蒜 散 量 五 7 本文 基于 地 面 气象 
观测 数据 ,利用 联合 国 粮农 组 织 (FAO ) 推 荐 的 Pen- 
man-Monteith 公式 计算 得 到 参考 作物 蒸 散 量 , 对 模 
拟 E7 进 行 长 时 间 的 尺度 扩展 ,计算 公式 如 下 : 

| 0408 A(R, - G) ey T u,(e,-e,) 


T,- " 
? A +y(1+0.34u,) (7) 


式 中 : ET, ABS EAR Be (mm d); A 为 饱和 
JK AUR Hl E BE (kPa 907) ; R, 为 地 表 净 辐射 (MJ， 
md ); C 为 土壤 热 通 量 (MJ.m…d 20; y 为 干 湿 
表 常 数 (kPa.% ;7 为 日 平均 气温 (%); u, 为 2m 
高 处 风速 (ms ; e, 为 饱和 水 气压 (kPa); e, HX 
际 水 气压 (kPa)。 
1.4.2 ARRE 结合 SEBAL 模 型 和 参考 作物 蒸 
散 量 计算 方法 对 月 实际 蒸 散 量 进 行 计算 : 

ET, 2 ET,F x ET, (8) 
XU ET, 为 实际 蒸 散 量 (mm) ; ETF 298 E ZR TX 
比 ,由 卫星 过 境 天 数 SEBAL 模型 计算 所 得 实际 蒸 散 
量 和 参考 燕 散 量 之 比 得 到 , 非 卫 星 过 境 天 数 的 
ET F 可 由 时 间 序 列 插值 方法 得 到 。 在 获得 连续 日 
实际 蒸 散 量 后 ,进行 累加 即 可 得 到 相应 月 蒸 散 量 。 
1.4.3 精度 评价 本 文 利 用 均 方 根 误差 (RHMSE ) 、 
决定 系数 ( 尼 ) 和 相对 误差 (MRE ) 等 定量 指标 对 模型 
反 演 进行 精度 验证 。 计 算 公 式 分 别 为 : 


(6) 


n 2 
RMSE- | y» (n. T ss] (9) 
1- 2. (E i -E S (10) 
EET) 


MRE = > (UE T -E T separ. iV E T d 


n 
式 中 ;n 为 观测 样本 个 数 ; ET, ,和 ETswar ;分 别 为 
用 于 验证 计算 所 得 ET 和 SEBAL 模 型 计算 所 得 7。 


NT (11) 
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1.5 水 量 平衡 及 灌溉 水 有 效 利用 系数 的 计算 

本 文通 过 对 玻 勒 河流 域 昌 马 灌 区 有 7 的 估算 , 结 
合 灌区 实际 灌溉 数据 和 降水 数据 ,进而 估算 得 到 洪 
区 在 不 同 灌溉 时 段 的 灌溉 水 有 效 利 用 系数 。 在 土 
壤 水 量 平衡 的 基础 上 ,规避 了 由 于 水 势 梯 度 作用 引 
起 的 水 分 渗 漏 和 补给 量 ,针对 灌区 建立 如 下 水 量 平 
衡 方程 ” : 

Aw z (P, * P, I) - (ET, +ET,+R) (12) 
式 中 : Aw 为 时 段 内 土壤 含水 量 的 变化 量 (m); P, 、 
P 忆 分别 为 时 段 内 非 灌溉 地 灌溉 地 降水 量 (m); 7 为 
时 段 内 的 净 灌 溉 量 (mm); ET, ET, 分 别 为 时 段 内 非 
灌溉 地 ,灌溉 地 蒸发 量 (mi); R 为 时 段 内 灌区 地 表 
径流 量 (m )。 

昌 马 灌区 地 处 干旱 的 平原 区 域 ,由 于 年 内 降水 
量 小 ,因此 忽略 降水 所 产生 的 地 表 径 流 ; 研 究 时 段 
内 土壤 层 仿 水 量变 化 较 小 ,本 研究 中 不 予 考虑 , 则 
可 对 灌区 内 的 水 量 平衡 方程 进行 简化 ,如 公式 
(13) 所 示 。 通 党 将 灌溉 地 消耗 的 灌溉 水 量 作为 洪 
溉 水 的 有 效 利 用 量 , 则 将 灌溉 水 有 效 利用 量 
(ET, - P,) 与 灌区 净 灌 溉 量 了 7 的 比值 定义 为 灌溉 水 有 
效 利 用 系数 7 , 见 公式 (14)。 

I- (ET, - P.) * (ET, - P.) (13) 

n=(ET,-P,)/1 (14) 

XB :CET,- P, ) 为 灌溉 地 消耗 的 灌溉 水 量 (m°); 

(ZE7,-P,) 非 灌溉 地 消耗 的 灌溉 水 量 (m^) ; m HE 
溉 水 有 效 利用 系数 。 


2 结果 分 析 


2.1 模型 验证 
本 文采 用 两 种 方法 验证 ET 模拟 结果 。 第 一 种 


T 


| (a) P-M 计 算 值 与 模拟 值 对 比 

6 = 
E 5 
"m 
m 
ES 3 
mM g RMSE = 0.49 mm * d! 

MRE = 8.46 % 
1 R? =0.87 
0 
0 2 4 6 8 
模拟 值 /mm 


是 采用 流域 内 4 个 气象 站 (敦煌 站 、 玉 门 站 、 马 坚 山 
站 、 瓜 州 站 ) 的 气象 观测 资料 ,利用 FAO P-M 公 式 计 
算得 到 各 站 点 4~10 月 的 日 参考 作物 蒸 散 量 ,结合 
生长 季 作 物 系 数 所 得 作物 实际 蒸 散 量 对 有 7 模拟 值 
进行 验证 ””。 第 二 种 方法 是 利用 蒸发 思 系 数 对 流 
域内 各 站 点 蒸发 亚 实 测 值 进行 折算 二 ,进而 与 了 7 模 
拟 结果 进行 对 比 验证 。 

上 述 两 种 方法 的 验证 结果 分 别 如 图 2a、2b 所 
IRo 在 提取 站 点 反 演 的 蒸 散 量 时 ,提取 以 站 点 为 中 
心 的 3 x 3 个 像 元 范围 内 的 平均 值 用 于 模型 验证 。 
从 图 2 可 知 ,两 种 验证 方法 的 尼 均 大 于 0. 8 ,表明 模 
型 反 演 值 与 P-M 计算 值 .实测 值 之 间 具 有 良好 的 相 
关 关 系 和 变化 趋势 ;两 种 方法 的 均 方 根 误差 相对 较 
小 , 且 0. 49 mm- d^! E H] P-M 计算 值 与 反 演 值 更 为 
接近 ;两 种 方法 的 相对 误差 均 为 8. 4% ,表明 模拟 值 
接近 实际 情况 ,模拟 结果 可 信 。 上 述 验 证 方法 的 结 
RRHH , SEBAL EH EAS A ROR Yt 49] 16] Pt LES) SE 
际 蒸 散发 情况 。 
2.2 蒸 散 发 模拟 结果 
2.2.1 日 蒜 散 量变 化 特征 图 3 分 别 为 模拟 得 到 
Hnc TAT Vitek 2017 4E 3 H ~ 20184E 2 AY H ZR ES 
间 分 布 ( 图 中 均 以 儒 略 日 DOY 表 示 日 期 ) ,(a) ~ (D 
分 别 代 表 日 序数 为 DOY009、DOY049、DOY089、 
DOY113, DOY129, DOYI61, DOY193, DOY225, 
DOY265 .DOY289 .DOY305 ,DOY337 HY H Z& #28 [i] 
分 布 。 

图 3 中 ,模拟 得 到 的 1 月 和 2 月 的 日 莹 散 量 低 ， 
分 别 为 0.58 mm*d ',1. 21 mm- d`’, CI EE RHE B 
量 以 低 值 为 主 ;3 ~5 月 的 日 蒸 散 量 逐 渐 增 大 ,分别 
达到 3. 00 mm-*d^',3. 52 mm:d^',4. 42 mm.d ,上 
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图 2 ET 模拟 结果 的 对 比 验 证 


Fig.2 Comparison and validation for simulated ET results 
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图 3 WAK A ZG B Al 
Fig.3 Distribution of daily evapotranspiration from 2017 to 2018 in the study area 


Tig DX Saas CERCHI] b ,高 值 区 即 红色 区 域 明显 增 
多 ;6 月 的 日 蒸 散 量 达 到 最 大 的 5. 03 mm*d …, 蒸 散 
量 高 值 区 主要 集中 于 上 游 地 区 , 且 中 下 游 增 幅 较 
大 ;7~10 月 的 日 均 蒸 散 量 开 始 下 降 , 但 7.8 月 的 日 
均 蒸 散 量 仍然 保持 一 个 较 高 的 蒸发 水 平 ,分 别 为 
4. 85 mm*d^',4. 62 mm:d ,其 中 9、10 月 的 日 均 蒸 
散 量 下 降 明 显 , 日 均 蒸 散 量 分 别 为 2. 70 mm'd ^, 
1. 62 mm.d ,其 中 上 游 地 区 蒸 散 高 值 区 降幅 较 大 ， 
区 域内 蒸 散 较 大 区 域 面积 明显 缩小 ;11~ 12 月 的 日 
蒸 散 量 到 达 年 内 最 低 值 ,分别 为 1. 1 mm*d^',0. 55 
mm:d ,区 域内 大 部 分 蒸 散 量 主要 分 布 在 0~1 
mm*“d-'。 总 体 来 说 ,研究 区 年 内 日 均 蒸 散 量 呈 先 增 
后 减 的 变化 趋势 ,最 大 值 出 现在 6 月 中 旬 ,最 小 值 出 
现在 12 月 下 旬 ,年 内 日 均 蒸 散 量 在 0. 55 ~ 5. 03 


mm: d` 内 变化 ,并 存在 明显 的 空间 分 布 差异 。 
2.2.2 PARE ”根据 对 2017 一 2018 年 的 
年 内 各 月 蒸 散 量 估算 ,获得 年 内 四 季 的 蒸 散 量 如 图 
4 所 示 , 图 中 (a) ~ (d) 分 别 代表 春季 (3 ~5 月 )、 夏 季 
(6~8 月 )、 秋 季 (9~11 月 )、 冬 季 (12~ 次 年 2 月 ) 四 
ZEE BU IH] AT ff DAD o 

MES 4 YF: ieee d) 18] o LAS [6] R5 2 ORE S 
显著 ,其 空间 分 布 特征 较为 相似 ,整体 表现 为 从 东 
南部 向 西北 部 逐渐 减少 。 春 季 蒸 散 量 空间 差异 较 
为 明显 ,上游 蒸 散 量 高 于 中 下 游 地 区 ,平均 燕 散 量 
为 165.51 mm。 这 是 由 于 春季 太阳 辐射 和 气温 的 逐 
渐 加 强 ,植被 在 生长 初期 的 蒸腾 作用 有 所 增强 。 进 
人 夏季 后 ,由 于 辐射 和 温度 达到 年 内 最 大 ,使 得 夏 
季 蒸 散 量 达到 年 内 最 大 的 201.83 mm, 从 而 使 得 蒸 
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散 量 的 空间 差异 最 为 显著 。 秋 季 气 温 \ 降 水量 和 太 
阳 辐 射 的 逐渐 减少 ,使 得 蒸 散 量 相对 夏季 呈现 明显 
的 下 降 , 流 域内 各 区 域 的 蒸 散 量 差 距 有 所 减 小 , 且 
中 下 游 地 区 的 葵 散 量 主要 处 于 200 mm 以 下 ,使 得 秋 
季 平 均 蒸 散 量 减 少 至 100. 51 mm, “FF IGA BOR 
为 四 季 内 最 小 , 仅 为 33. 92 mm, 冬季 蒸 散 量 的 空间 
差异 最 小 , 仅 上 游 部 分 山区 蒙 散 量 稍 高 ,其 余 区 域 
均 为 小 于 50 mm 的 低 值 区 即 蓝 色 部 分 。 因 此 ,年 内 
季节 蒜 散 量变 化 为 :夏季 > 春季 > 秋季 > 冬季 。 
结合 图 5 的 研究 区 土地 利用 状况 可 以 看 出 ,四 
季 的 蒸 散 量 高 值 区 域 主 要 集中 在 上 游 地 区 ,分 布 在 
紧邻 山区 的 草地 林地、 河流 沿岸 和 部 分 水 体 区 域 ; 
中 下 游 地 区 燕 散 量 处 于 整个 流域 的 较 低 水 平 , 蒸 散 


ET/ mm 0 


主要 集中 于 中 下 游 地 区 中 的 低 覆 盖 度 草地 , 蒸 散 量 
较 小 ;流域 西部 地 区 是 整个 流域 蒸 散 量 最 小 的 区 
域 , 主 要 与 中 游 地 区 大 面积 分 布 的 沙 地 .戈壁 和 未 
利用 土地 相关 ,具有 极 低 的 蒸 散 量 , 则 人 研究 区 蒸 散 
量 大 小 在 空间 上 呈现 上 游 明 显 大 于 中 下 游 地 区 。 
2.3 水量 平衡 及 灌溉 水 有 效 利用 系数 计算 结果 

本 文 对 玻 勒 河流 域内 面积 最 大 的 昌 马 灌区 采 
用 遥感 蒸 散 发 模型 ,进行 灌溉 水 有 效 利 用 系数 估算 
与 灌溉 效率 评价 。 昌 马 灌区 位 于 流域 中 游 地 区 , E 
要 分 布 有 耕地 和 部 分 林地 (图 1) ,其 年 内 灌溉 时 段 
划分 为 :春灌 (3 ~4 月 ) 、 夏 灌 (5~7 月 ) 秋江 (8~9 
月 ) 冬 灌 (10~11 月 )。 昌 马 灌区 春灌 BWE PATE 
和 冬 灌 的 ET 结果 如 图 5 中 (a) ~ (d) 所 示 。 灌 区 中 


50 100 150 200 250 300 350 400 
图 4 研究 区 四 季 蒸 散 量 空间 分 布 图 


Fig.4 Distribution of seasonal evapotranspiration from 2017 to 2018 in the study area 
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5 昌 马 灌区 各 灌 季 蒸 散 量 分 布 图 


Fig.5 Distribution of seasonal evapotranspiration in Changma irrigated area 
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由 于 各 灌 季 的 灌溉 用 水 量 和 种 植 情况 的 不 同 ,使 得 
各 灌 季 的 蒸 散 量 存在 明显 的 时 空 差 异 ,其 中 5~7 月 
的 夏 灌 阶段 #7 最 大 ,其 次 为 秋 灌 和 春灌 阶段 , 冬 灌 
阶段 的 ET 最 小 ;从 ET 空间 分 布 可 以 看 出 ,不同 灌溉 
时 段 内 的 £7 高 值 区 主要 分 布 在 中 部 与 东南 部 , 低 值 
区 主要 分 布 在 边缘 区 域 。 

结合 实际 灌区 内 降水 数据 .灌溉 面积 和 净 灌溉 
量 数据 计算 水 量 平衡 所 需 分 量 , 其 中 灌区 3 ~ 11 月 
ET 及 降水 量变 化 如 图 6 所 示 。ET 变 化 呈 先 增 后 减 
的 变化 趋势 ,3 ~ 6 月 逐渐 增 大 ,并 在 6 月 达到 最 大 的 
129 mm;6~11 月 逐渐 下 降 , 并 在 11 月 达到 最 小 的 
20 mm。3~11 月 的 各 月 降水 量变 化 幅度 较 大 ,其 中 
6.7 月 降水 量 最 高 , 均 大 于 30 mm;4 月 和 8 月 的 降水 
量 较 大 , 均 在 15 mm 左右 ;其 余 月 份 的 降水 量 均 处 于 
较 低 水 平 , 均 小 于 7 mm。 上 述 的 各 月 蒸 散 量 和 降水 
量 作为 水 量 平衡 的 重要 分 量 , 将 对 不 同 灌 季 的 灌溉 
效率 产生 较 大 的 影响 。 

通过 对 灌溉 水 有 效 利 用 量 的 计算 ,结合 各 灌 季 
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Fig.6 Variation of monthly evapotranspiration and 


precipitation in Changma irrigated area 
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图 7 昌 马 灌区 各 灌 季 参数 变化 


Fig.7 Variation of monthly parameters in Changma 


irrigated area 


的 灌溉 面积 和 净 灌 溉 量 , 求 得 灌区 内 各 灌 季 的 灌溉 
水 有 效 利用 系数 ,结果 如 图 7 所 示 。 计 算 所 得 的 春 
WE RE 、 秋 灌 和 冬 灌 的 灌溉 水 有 效 利 用 系数 分 别 
为 0. 76,0. 71,0. 69,0. 55 ,年 内 平均 灌溉 水 有 效 利 
用 系数 为 0. 67 ,其 中 冬 灌 的 灌溉 水 有 效 利 用 系数 低 
于 年 平均 水 平 ; 灌 溉 水 有 效 利用 系数 在 年 内 旦 逐渐 
减少 的 趋势 ,其 中 春灌 的 灌溉 水 有 效 利用 系数 最 
大 ,其 次 为 夏 灌 和 秋 灌 , 冬 灌 最 低 。 年 内 净 灌 溉 量 
和 灌溉 水 有 效 利 用 量 差 异 显著 ,其 中 春灌 用 水 量 最 
小 , 夏 灌 最 大 ,其 次 为 秋 灌 和 冬 灌 ,这 与 年 内 不 同 季 
节 的 作物 种 植 面积 及 作物 生长 阶段 的 需 水 量 关系 


3 讨论 


3.1 BARRED RIB 

本 文 利用 MODIS ERAAI SEBAL TUS! , OTL 
勒 河流 域 2017 一 2018 4F. PY B ZS Rich BE TT E Rt fh SA- 
并 分 析 其 空间 分 布 特征 。 结 果 表 明 , 琉 勒 河流 域 年 
内 日 均 蒸 散 量 呈 单 峰 变 化 趋势 ,并 在 0. 55 ~5. 03 
mm*d 内 变化 ,四 季 中 夏季 的 蒸 散 量 最 大 ,春秋 次 
之 ,冬季 最 小 。CHANG 5E" H] SEBAL 模型 分 别 
5 i8 C SL Egr6HmHzutig4r T 3-5 
mm: d ', JA EJES 7 38 et SEBAL E70 Fz 18615: SU tt 
勒 河流 域 8 月 多 年 平均 日 蒸 散 量 为 0~ 8. 52mm: 
d-', 上 述 结论 与 本 文 所 得 相应 时 段 结果 相近 ,证 明 
本 文 计 算 结 果 人 合理。 年 内 莹 散 量 的 变化 由 于 气温 、 
降水 ,太阳 辐射 .植被 生长 态势 和 土壤 含水 量 等 因 
素 的 共同 影响 ,并 表现 出 季节 性 变化 。 如 夏季 时 段 
是 植被 生长 旺季 阶段 ,充足 的 水 分 和 强 辐射 使 得 水 
体 和 植被 的 蒸发 和 蒸腾 作用 达到 最 强 ,而 冬季 蒸 散 
由 于 气温 和 辐射 低 ,降水 少 , 地 表 植 被 和 土壤 水 分 
蒸发 量 较 小 ,造成 冬季 蒸 散 量 低 。 

从 反 演 结果 上 看 , 莹 散发 的 空间 分 布 与 土地 履 
被 类 型 具有 高 度 一 致 性 。 蒸 散发 高 值 区 主要 分 布 
在 琉 勒 河上 游 地 区 高 履 盖 度 草 地 以 及 水 域 周围 等 
区 域 ,这 是 由 于 上 游 地 区 草地 有 来 自 于 山区 季节 性 
降水 的 水 源 补充 ,土壤 湿润 ,利于 蔡 发 。 玻 勒 河流 
域 中 下 游 的 沙 地 与 荒漠 分 布 广泛 ,并 有 小 面积 绿洲 
交错 分 布 ,为 流域 BT7 低 值 的 主要 分 布 区 域 。 这 是 由 
于 荡 漠 区 草地 和 林地 的 地 下 水 资源 极其 稀缺 ,补给 
水 源 主要 依赖 于 地 下 水 和 降水 的 补给 ,加 之 沙 地 、 
戈壁 区 域 水 资源 匮乏 .地表 裸露 ,共同 导致 了 该 区 
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域 的 蒸 散 量 低 。 这 与 金 学 杰 等 光 对 干旱 区 黑河 流 
域 的 蒸 散发 估算 结果 中 水 体 和 草地 蒸 散 发 量 大 ,而 
沙 地 和 裸 地 的 蒸 散 量 最 小 的 结论 一 致 。 
3.2 灌溉 水 有 效 利 用 系数 计算 

AS SCF FE FE ERR OY OP PTT Yat Sak HE 
区 的 灌溉 水 有 效 利 用 系数 进行 了 估算 ,用 于 计算 的 
降水 量 .实际 灌溉 面积 和 净 灌 溉 量 均 为 获取 的 实测 
统计 数据 ,因此 所 得 灌溉 水 有 效 利用 系数 具有 一 定 
的 精度 。 本 文 计算 得 到 2017 年 平均 灌溉 水 有 效 利 
用 系数 为 0. 67 ,不 同 灌溉 时 段 的 灌溉 水 有 效 利 用 系 
数 存 在 一 定 差异 。 春 灌 由 于 处 于 作物 生长 初期 ,其 
灌溉 用 水 量 较 小 , 莱 发 量 较 大 ,使 得 灌溉 水 有 效 利 
用 系数 较 大 ; 夏 灌 和 秋 灌 的 灌溉 水 有 效 利 用 系数 相 
对 减少 ,这 是 由 于 夏季 和 秋季 阶段 为 满足 作物 生长 
旺盛 期 的 需 水 量 , 净 灌 溉 量 达 到 最 大 ,但 过 度 灌 溉 
会 导致 灌溉 水 在 运输 过 程 中 经 过 渠道 及 田间 的 水 
量 损失 和 渗 漏 , 且 该 时 段 的 降水 量 较 大 ,从 而 使 得 
该 阶段 的 灌溉 水 有 效 利 用 系数 下 降 ,因此 可 通过 减 
少 灌溉 过 程 中 水 量 输送 的 损失 来 提高 灌溉 效率 ; 冬 
灌 的 灌溉 水 有 效 利 用 系数 相对 较 低 ,主要 是 因为 灌 
区 内 10~ 11 月 处 于 作物 生长 未 期 ,灌溉 地 的 蒸 散 量 
迅速 减少 ,使 得 该 阶段 灌溉 水 有 效 利 用 系数 有 所 
下 降 。 

灌溉 水 有 效 利用 系数 同时 受到 净 灌 溉 水 量 、 降 
水 量 和 ET 的 影响 ,也 与 灌溉 技术 和 输 水 系统 相关 。 
因此 可 根据 不 同 季 节 的 蔡 散 特征 和 灌溉 水 利用 情 
况 , 采 取 不 同 的 节 水 灌溉 措施 和 合理 的 灌溉 水 量 来 
提高 作物 的 灌溉 效率 。ET 作 为 灌溉 耗 水 和 灌溉 水 
有 效 利用 系数 计算 的 重要 组 成 ,E7 的 准确 估算 对 认 
识 作物 蒸 散 特 征 和 耗 水 情况 具有 重要 影响 。 因 此 ， 
在 今后 的 研究 中 应 根据 实际 情况 改善 模型 参数 , 同 
时 结合 高 分 辩 率 影像 将 有 效 提高 流域 及 灌区 MUS] 
反 演 精度 ,进而 更 加 准确 估算 灌区 的 灌溉 水 利用 
情况 。 


4 结论 


v 


(1) ifii By 19] jfi ER 2017—2018 年 年 内 平均 日 
蒸 散 量 总 体 呈 单 峰 变化 趋势 ,年 内 日 均 蒸 散 量 
在 0.55~5.03 mm-d 内 波动 变化 。 日 燕 散 量 从 3 
月 开始 波动 上 升 ,并 在 6 月 达到 最 大 值 5. 03 mm 
d ,之 后 的 7 月 ~ 次 年 的 2 月 的 日 落 散 量 呈 波动 下 


降 趋 势 , 其 中 日 蒸 散 量 最 小 值 0. 55 mm .d 出 现在 
12 月 至 次 年 1 月 之 间 。 

(2) iene 30] Ye Je DUE 7 变化 显著 ,四 季 Z7 值 
从 大 到 小 排列 为 :夏季 (201. 83 mm) > 春季 (165. 51 
mm) > 秋季 (100. 51 mm) > 冬季 (53. 92 mm), 。 四 
Æ 7 的 时 空 分 布 异 质 性 显著 ,时 间 上 ,夏季 和 秋季 
的 57 差异 较 春 季 和 冬季 明显 ;空间 上 aito 9) do] D RK 
ET 由 东南 向 西北 逐渐 减 小 ,上 游 地 区 ET 明显 高 于 
中 下 游 地 区 , 藻 散 发 的 空间 分 布 与 土地 覆 被 类 型 具 
有 高 度 一 致 性 。 

(3) 玻 勒 河流 域 昌 马 灌区 由 于 各 灌 季 的 灌溉 用 
水 量 和 种 植 情况 的 不 同 ,不 同 灌溉 时 段 的 灌溉 地 薄 
散 量 存在 明显 的 空间 分 布 差 异 与 空间 变化 ,灌区 内 
KEER ET 高 值 区 主要 分 布 在 中 部 与 东南 部 , 低 
值 区 主要 分 布 在 西北 部 和 灌区 边缘 。 

(4) 根据 水 量 平衡 对 昌 马 灌区 灌溉 水 有 效 利用 
系数 进行 估算 ,降水 量 与 蒸 散 量 在 治 季 内 呈 先 增 后 
减 的 波动 变化 。 昌 马 灌 区 年 内 灌溉 水 有 效 利 用 系 
数 呈 下 降 趋势 ,计算 所 得 的 春灌 、 夏 灌 、 秋 灌 和 冬 灌 
各 灌 季 的 灌溉 水 有 效 利用 系数 分 别 为 0. 76 .0. 71. 
0. 69 .0. 55 ,年 内 平均 灌溉 水 有 效 利 用 系数 为 0. 67。 
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Estimation of evapotranspiration in Shule River Basin based on 
SEBAL model and evaluation on irrigation efficiency 


NING Ya-zhou', ZHANG Fu-ping, FENG Qi, WEI Yong-fen’, 
LI Ling', LIU Jie-yao', ZENG Pan-ru' 
(1 School of Geography and Tourism ,Shaanxi Normal University , Xi' an 710119 , Shaanxi, China; 


2 Northwest Institute of Eco-Environment and Resources , Chinese Academy of Sciences , Lanzhou 730000, Gansu , China; 


3 River Basin Research Center , Gifu University , Gifu 501-1193 , Japan ) 


Abstract: Water resources are scarce for inland river basins in arid regions of northwestern China, where water 
consumption is mainly for agriculture. Especially in arid regions , evapotranspiration ET) is the main link to the wa- 
ter cycle and has an important impact on environmental changes. It is therefore important to accurately estimate ET 
and irrigation efficiency for predicting climate change and for rationally allocating water resources in inland river ba- 
sins. The Shule River Basin is one of three inland river basins in the Hexi Corridor, and its ecological environment 
is fragile given its large annual average ET. Water use accounts for a large proportion , and the utilization rate of ag- 
ricultural irrigation water is high in this region. In this study, we estimated the daily evapotranspiration (ET) of 
the Shule River Basin (SRB) based on the Surface Energy Balance Algorithms for Land (SEBAL) model with 
satellite data (from the moderate-resolution imaging spectroradiometer known as MODIS) and meteorological 
data from March in 2017 through February in 2018. Then we applied a time-scale extension to obtain monthly 
and seasonal ET by using daily ET and meteorological data. We also analyzed the spatiotemporal distribution 
characteristics of ET to derive daily and seasonal ET. Then, combining the ET results from the SEBAL simulation 
and data collected on irrigation and precipitation for the Changma irrigation district (CM) during 2017 — 2018, 
we calculated the coefficient of irrigation-water effective utilization (“n” ) for the CM on the basis of the wa- 
ter-balance equation, which was used for irrigation-efficiency assessment for the CM. 

The main conclusions are as follows: (1) The variation range of average daily ET values was 0.55 — 5.03 
mm across the SRB during 2017 — 2018, which showed a single peak change and obvious difference in spatial dis- 
tribution. (2) There were significant differences in seasonal ET in the SRB ; summer ET was the highest (201.83 
mm) , followed by spring and autumn, and the winter ET was the lowest (53.92 mm). The spatial pattern of sea- 
sonal ET shows that ET in the upper reaches of SRB was significantly higher than that in the middle and lower 
reaches , and ET decreased from the southeast to the northwest in the SRB. (3) The quantity of irrigation water in 
different irrigation periods resulted in the ET differences in the CM. The high-ET area was mainly in the middle 
and southeastern parts of the irrigation district, and the low-ET area was mainly in the northwestern part and 
along the margins of the district. (4) The “n” of CM showed a downward trend during the irrigation period, in 
which the values of “n” for the spring, summer, autumn, and winter irrigation periods were 0.76 ,0.71, 0.69, and 
0.55 , respectively. The CM had a higher irrigation efficiency , with the annual “n” of 0.67. 

Through simulation of monthly ET in the SRB for the period 2017 — 2018, the annual ET spatial distribu- 
tions and fluctuations were more detailed and clearer. We found that ET of irrigation land could be accurately esti- 
mated by remote-sensing and meteorological data. The “n” of the irrigation district could easily be obtained from 
ET and measured irrigation data in each irrigation period, which can reflect CM irrigation efficiency during the 
year in the SRB. 

Key words:  evapotranspiration; SEBAL model; Shule River Basin; coefficient of irrigation water effective 
utilization; Changma irrigation district 


